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ABSTRACT 

We exploit the galaxy groups catalogue of Yang et al. and the galaxy properties mea- 
sured in the SDSS Data Releases 4 and 7 to study how the gas-phase metallicities of 
Q^) ' star- forming galaxies depend on environment. We find that satellite and central galax- 

ies follow a qualitatively similar stellar mass (M+) - gas-phase metallicity relation, 
whereby their gas-phase metallicity increases with M t . Satellites, though, have higher 
I | gas-phase metallicities than equally massive centrals, and this difference increases with 

decreasing stellar mass. We find a maximum offset of 0.06 dex at log(M i , / h~ 2 M Q ) ~ 
8.25. At fixed halo mass, centrals are more metal rich than satellites by ~0.5 dex 
on average. This is simply due to the fact that, by definition, centrals are the most 
massive galaxies in their groups, and the fact that gas-phase metallicity increases 
with stellar mass. More interestingly, we also find that the gas-phase metallicity of 
satellites increases with halo mass (Mh) at fixed stellar mass. This increment is more 
pronounced for less massive galaxies, and, at M* ~ 1O 9 /i _2 M , corresponds to ~ 0.15 
dex across the range 11 < log(Mh//i _1 M Q ) < 14. We also show that low mass satellite 
galaxies have higher gas-phase metallicities than central galaxies of the same stellar 
metallicity. This difference becomes negligible for more massive galaxies of roughly 
solar metallicity. We demonstrate that the observed differences in gas-phase metal- 
^^O . licity between centrals and satellites at fixed M+ are not a consequence of stellar 

mass stripping (advocated by Pasquali et al. in order to explain similar differences 
but in stellar metallicity) , nor to the past star formation history of these galaxies as 
quantified by their surface mass density or gas mass fraction. Rather, we argue that 
these trends probably originate from a combination of three environmental effects: (i) 
strangulation, which prevents satellite galaxies from accreting new, low metallicity gas 
which would otherwise dilute their ISM, (ii) ram-pressure stripping of the outer gas 
| disk, thereby inhibiting radial inflows of low-metallicity gas, and (iii) external pres- 

sure provided by the hot gas of the host halo which prevents metal-enriched outflows 
from escaping the galaxies. Each of these three mechanisms naturally explains why 
the difference in gas-phase metallicity between centrals and satellites increases with 
decreasing stellar mass and with increasing host halo mass, at least qualitatively. How- 
ever, more detailed simulations and observations are required in order to discriminate 
between these mechanisms, and to test, in detail, whether they are consistent with 
the data. 

Key words: galaxies: abundances - galaxies: evolution - galaxies: fundamental pa- 
rameters - galaxies: groups: general - galaxies: star formation 
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1 INTRODUCTION 

Star forming galaxies in the local Universe define a funda- 
mental plane, where their star formation rate (SFR) and 
gas-phase metallicity (typically 12 + log(0/H)) correlate 
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with their stellar mass (M*, cf. Lara-Lopez et al. 2010; Man- 
nucci et al. 2010). The individual r elations in this plane 
have been long known; for example, iLequeux et all (|l979l ) 
were the first to recognize the existence of a correlation be- 
tween galaxy magnitude and gas-phase metallicity, whereby 
more luminous galaxies exhibit higher metallicities. More 
observational evidence in support of such a dependence was 
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later provided by Garnett & Shields (1987), Skillman et 
al. (1989) and Zaritsky et al. (1994) among others. The 
advent of t he Sloan Digita l Sky Survey (SDSS) made it 
possible for iTremonti et all (|2004) to transform the galaxy 
luminosity - metallicity relation into the more fundamen- 
tal dependence of gas-phase metallicity on galaxy stellar 
mass, where the interstellar medium (IS M) of more massive 
galaxi es is metal richer. Using the SDSS. lBrinchmann et all 
(2004) and Salim et al.l (|2005T ) established a correlation be- 
tween the SFR of galaxies and their stellar mass, whereby 
more massive, star forming galaxies sustain larger star for- 
mation rates (see also Ellison et al. 2008) . The star formation 
rate has been seen to correlate also with the gas-phase metal- 
licity of galaxies, although in a way that dif fers according to 
the ga laxy stellar mass. As pointed out bv lMannucci et ahl 
|2010l ). the gas-phase metallicity of low-mass galaxies de- 
creases as their SFR increases, while high-mass galaxies 
do not show any significant dependence of their gas-phase 
metallicity on SFR (see also Yates et al. 2012). 

Such an interplay among Af* , SFR and gas-phase metal- 
licity of galaxies is usually explained in terms of stellar feed- 
back triggering gas and metals outflows (see Edmunds 1990, 
Lehnert & Heckman 1996; Frye, Broadhurst & Bemtez 2002; 
Garnett 2002; Tremonti et al. 2004; Kobayashi, Springel 
& White 2007; Scannapieco et al. 2008; Weiner et al. 
2009; Spitoni et al. 2010; McCarthy et al. 2011; Peeples & 
Shankar2011).|M annucci et al. I (l20ia ) envisaged a condition 
of steady-state in the local Universe, where infall of metal- 
poor gas (diluting the galaxy gas-phase metallicity and sus- 
taining its star formation activity) occurs together with out- 
flow of metal-rich gas whose efficiency should depend on the 
galaxy mass and SFR. Alternative explanations of the fun- 
damental plane of star-forming galaxies invoke downsizing 
(where the star formation efficiency depends on galaxy mass, 
e.g. Brooks et al. 2007; Mouhcine et al. 2008; Calura et al. 

2009) , a dependence of the initial mass function on galaxy 
mass (Koppen, Weidner & Kroupa 2007) or a model of pure 
infall of metal-poor gas (Finlator & Dave 2008; Dave et al. 

2010) . 

One can easily expect that environment may intervene 
in shaping the fundamental plane of star-forming galaxies. 
Environmental processes, such as strangulation (i.e. the re- 
moval of the gas reservoir of satellites after they are ac- 
creted by their host halo), ram pressure stripping (i.e. the 
removal of gas from satellites moving through a dense intra- 
cluster medium), galaxy harassement (i.e. high-speed impul- 
sive galaxy encounters) and galactic wind confinement (i.e. 
due to the pressure exerted by the hot gaseous atmosphere of 
the host halo, which can inhibit galactic winds) , are believed 
to alter the inflows/outflows of gas experienced by galaxies, 
modifying their star formation activity and their gas-phase 
metallicities with respect to that predicted by closed-box 
evolution. 

Observations of the ISM in galaxies residing in the 
Lynx-Cancer void show that the gas-phase metallicity of 
these galaxies is lower by 30% on average than that 
of equally-massive galaxies in higher-density environments 
(Pustilnik, Tepliakova & Kniazev 2011). 

Nearby clusters (such as Coma, Hercules and Virgo) 
have been extensively observed in search of environ- 
mental effects o n the properties of star-forming g alax- 
ies. For example, IShields. Skillman fc Kennicuttl |l99l[ ) and 



ISkillman et al] (|l996h found that the ISM of spiral galax- 
ies in Virgo is on average metal-richer than that of com- 
parable field galaxies. When a distinction is made on the 
basis of their HI content, the Virgo spirals that are HI de- 
ficient and also located closer to the cluster core are seen 
to have a higher gas-phase metallicity and a significantly 
lower birthrate (i.e. the ratio of newly-born stars to the stars 
formed in the past; Gavazzi et al. 2002) than those with 
normal HI content, which are usually located in the clus- 
ter outskirts a nd have similar me t allicit y as field galaxies. 
More recently, IPetropoulou et all (|201ll ) have shown that 
spiral galaxies in the Hercules cluster are chemically evolved 
with an oxygen abundance close to solar, and their gas-phase 
metallicity does not depend on the local density within the 
cluster. Interestingly, their Ha emission is less spatially ex- 
tended than their optical disks and their star formation ac- 
tivity takes place preferentially in their inner regions. This 
result supports the picture where ram-pressure stripping due 
to the intergalactic medium (IGM) has truncated the HI 
disk of these spirals, with the result of preventing inflows of 
metal-poor gas from the outer disk and halo which would 
otherwise dilute the gas-phase metallicity in the inner re- 
gions of these galaxies. Suppression of infall of metal-poor 
gas due to environment would also be an explanation for the 
higher gas-phase metallicity of cluster spirals. 

As for dwarf galaxies, th e studies of IVflche3 (| 19951 ) and 
iLee. McCall fc Richerl (|2003l ) did not reveal a significant de- 
pendence of gas-phase metallicity on environment. On aver- 
age, equally-bright dwarfs residing in the field and in Virgo 
do not differ in their O/H abundance, but some Virgo dwarfs 
can show a gas deficiency when compared with equally 
metal-rich dwarfs in the field. Such a gas deficiency seems to 
correlate with the X-ray surface brightness of the IGM, thus 
suggesting that ram-pressure stripping due to the IGM de- 
prived these Virgo dwarfs of some of their gas. Ram-pressure 
stripping has been observed 'in action' for already a number 
of galaxies in Virgo (e.g. Crowl et al. 2010), whose HI disks 
appear distorted and truncated with respect to their stel- 
lar disks. A somew hat stronger dependence o n environment 
was established by IPetropoulou et ahl |201ll ) for the dwarf 
galaxies in the Hercules cluster, whose gas-phase metallicity 
is seen increasing with local density. Metal-poorer dwarfs are 
associated with lower densities (typically the infall regions of 
the cluster) and with a stronger Ha emission with respect to 
the global dwarf population in the cluster (Mahajan, Haines 
& Raychaudhury 2010), a sign that they undergo a starburst 
phase upon accretion onto the cluster. Such a star formation 
activity has to get quenched at later times in order to repro- 
duce the observed fractions of red-and-dead satellites (e.g. 
van den Bosch et al. 2008; Kang & van den Bosch 2008; 
Kimm et al. 2009; Pasquali et al. 2009; Weinmann et al. 
2009). 

The large statistical power of the SDSS survey has 
not yet provided a definitive assessment of the envi- 
ronmental dependence of galaxy ga s-phase metallicity. In 
iMouhcine. Baldrv fc Bamfordl (|2007l ) the gas-phase metal- 
licity - mass relation turns out to be weakly dependent 
on environment, defined as the number density of galax- 
ies within the projected distance to the closest 4th and 
5th neighbours. Star-forming galaxies show a marginal in- 
crease of their O/H abundance at fixed stellar mass in 
denser environments, equivalent to a few percent for the 
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Table 1. Samples statistics 







Description 




Centrals 


Satellites 


Satellite Fraction 


Sample 


S 


galaxies with spectral S/N ^ 20 and measured af 


;e and Z 


70,067 


13,626 


19.4% 


Sample 


G 


galaxies with EW (H/3) ^ -3 A and measured 12 


+ log(0/H) 


69,688 


14,182 


20.4% 


Sample 


C 


galaxies that are both in Samples S and G 




8,134 


1,595 


19.6% 



more massive galaxies up to ~20% at log(M*/Mo) < 
9.5. The authors concluded that the evolution of star 
forming galaxies is mainly driven by their intrisic prop- 
erties, and does n ot depend on environment at large. 
ICooper et al.l (2008) found a stronger dependence of gas- 
phase metallicity on galaxy environment (denned as the 
number density of galaxies within the projected distance 
to the closest 3rd neighbour), in that galaxies residing in 
higher density environments are metal-richer (see also El- 
lison et al. 2009 ). The question arises whethe r these dif- 
ferences betw een Mouhcine. Baldrv fe Bamfordl (|2007h and 
ICooper et al.l (2008?) may be due to differ ent definitions of 
enviro nment. As alre ady discussed b y e.g. [ Weinmann et all 
l|200rj ) and recently iMuldrew et all 120111 ). the projected 
number density of galaxies is of 'tricky' interpretation as 
it depends on environment itself. In massive haloes such as 
clusters, it probes environment on a scale smaller than the 
halo virial radius, while in low mass haloes the environment 
probed typically encompasses multiple dark matter haloes. 



In this paper, we revisit the issue of gas-phase metallic- 
ity vs. galaxy environment by making use of the catalogue 
of galaxy groups ext racted from the SDSS Data Release 4 
bv lYang et ail (|2007f) . In this catalogue, environment is de- 
fined as the total amount of dark matter associated with 
each galaxy group, and also allows us to distinguish (and 
treat separately) central from satellite galaxies, for which 
semi-analytic models of galaxy formation and evolution pre- 
dict rather different evolutionary paths. We aim at achiev- 
ing a detailed assessment of how gas-phase metallicities of 
galaxies depend on galaxy stellar mass (over the interval 8 
< log(M*//i -2 M0) < 11), hierarchy (centrals vs satellites) 
and environment across a large range of halo masses (11 < 
log(Mh//i _ Mo) < 14), in order to gain important insights 
into how environment regulates gas inflows feeding the star 
formation activity of galaxies, and affects the efficiency of 
galactic outflows through the pressure applied by the intra- 
group hot gas. 



This paper is organized as follows. In Sect. 2 we describe 
the c atalogue of galaxy groups obtained by lYang et al.l 
(120071 1 for the SDSS Data Release 4 (DR4). The stellar mass 
- gas-phase metallicity relation obtained for satellites and 
centrals separately is discussed in Sect. 3, while the halo 
mass - gas-phase metallicity relation is presented in Sect. 4. 
The results are discussed in Sect. 5 and the conclusions fol- 
low in Sect. 6. Throughout this paper we adopt a flat ACDM 
cosmology with fi m = 0.238 and = 0.762 (Spergel et al. 
2007) and we express units that depend on the Hubble con- 
stant in terms of h = _r7 o /(100 km s _1 Mpc -1 ). 



2 DATA 

Our analysis is based o n the SDSS DR4 c atalogue of galaxy 
groups constructed by lYang et al.l (12007 ) by applying the 
halo-based group finder algorithm of lYang et al.l (|2005T ) to 
the New York University Value-Added Galaxy Catalogue 
(NYU-VAGC, Blanton et al. 2005) extracted from SDSS 
DR4 (Adelman-McCarthy e t al. 2006) . From the NYU- 
VAGC Main Galaxy Sample lYang et al] (|2007l l selected all 
galaxies with an apparent magnitude (corrected for Galac- 
tic foreground extinction) brighter than r = 18 mag, in the 
redshift range 0.01 ^ z ^ 0.20 and with a redshift com- 
pleteness C z > 0.7. These galaxies were used to build three 
galaxy groups samples: sample I, which only uses the 362356 
galaxies with measured redshifts from the SDSS; sample II, 
which includes an additional 7091 galaxies with SDSS pho- 
tometry but redshifts taken from alternative surveys; sample 
III which lists an additional 38672 galaxies lacking a redshift 
due to fibre collisions, but being assigned the redshift of their 
nearest neighbour (cf. Zehavi et al. 2002). This paper focuses 
solely on sample II, where galaxies are split into 'centrals' 
(defined as the most massive group members on the basis of 
their M*) and 'satellites' (all group members that are not 
centrals). Galaxy magnitudes and colours are based on the 
standard SDSS Petrosian technique (Petrosian 1976; Strauss 
et al. 2002), are corrected for Galactic foreground extinction 
(Schlegel, Finkbeiner & Davis 1998) and J^-corrected and 
evolution corrected t o z = 0.1 according to the procedure of 
iBlanton et ail (|2003h . We use the notation °' 1 Mx to indi- 
cate the resulting absolute magnitude in the photometric X 
band. Galaxy stellar masses are computed using t he relation 
betwe en stellar mass-to-light ratio and colour of iBell et al.l 
(2003), and galaxy morphology is defined in terms of the 
concentration parameter (C — rgo/rso where rgo and rso 
are the radii encompassing 90% and 50% of the Petrosian 
r— band flux, respectively). Values of C larger than 2.6 typ- 
ically refer to bulge-dominated galaxies (cf. Strateva et al. 
2001). 

For each group in the catalogue I Yang et alj (|2007l ) pro- 
vided two estimates of its dark matter halo mass, Mh, one 
based on the ranking of its total luminosity, the other on 
the ranking of its total stellar mass. These halo masses de- 
viate from each other by about 0.1 dex at l ow Mh and 
0.05 d ex at the massive end. The method of I Yang et al.l 
(2007) is able to assign Mh only to groups more massive 
than ~ 10 12 /i _1 Mq and with one or more members brighter 
than 01 M r — 51ogft = -19.5 mag. For all the other groups, 
lYang. Mo fc van den Boschl l|2008l ) used the relations be- 
tween the luminosity (or M*) of central galaxies and the 
halo mass of their groups to estimate the halo mass of single 
central galaxies down to Mh ~ 10 1:L /i _1 Mq. In what fol- 
lows we will use the halo masses derive d from the group 's 
total stellar mass, because, as shown bv lMore et ail l|201lft , 
stellar mass is a better indicator of halo mass than luminos- 
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Figure 1. The normalized distributions of sample II (black dot-dashed line) and samples S and G (solid grey and black, respectively) 
in colour, concentration, stellar mass and halo mass. The distributions are not weighted by 1/V m ax- 
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Figure 2. The normalized distributions of centrals (grey) and satellites (black) in sample G in stellar mass, gas-phase metallicity and 
halo mass. The distributions are not weighted by 1/Knax- 



ity. Since the groups catalogue of I Yang et alj |2007l ) is not 
volume limited, it suffers from the Malmquist bias, which 
causes an artificial increase of the average luminosity (and 
M^) of galaxies with increasing redshift. This effect is espe- 
cially important for satellites, since they span a large range 
of stellar masses within each halo. We take the Malmquist 
bias into account by weighting each galaxy by 1/V ma x, where 
V m ax is the comoving volume of the Universe out to a co- 
moving distance at which a galaxy would still have made the 
selection criteria of the groups sample. Therefore, all distri- 
butions in the following sections are weighted by 1/V max 
unless specified otherwise. 



2.1 Stellar population parameters: sample S 

We matched sample II w ith the catalogue of s tellar ages and 
metallicities obtained bv lGallazzi et al] ll2005l) f or the galax- 
ies in the SDSS DR4. Briefly, iGallazzi et al.l l|2005h com- 
pared the strength of stellar absorptions in the observed 
spectrum of each galaxy with the predictions of a Monte 
Carlo library of 150 000 star formatio n histories (SFHs), 
based on the lBruzual fc Charlotl (|2003l) population synthe- 
sis code and the IChabrierl 1 20031 ) initial mass function. The 
SFHs in the library have a declining star formation rate 
(SFR) over varying timescales, and are superposed with ran- 
dom bursts of varying intensity and duration. These bursts 
are added in such a way that only 10% of the models can ex- 
perience a burst of star formation in the last 2 Gyr. For each 
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galaxy a probability density distribution of stellar r— band 
flux-weighted stellar ages and metallicities can be derived, 
whose median values we will use for our analysis below. It 
has to be noted that the stellar ages and metallicities de- 
rived as above refer to the redshift at which galaxies are 
observed. No correction to z — was attempted because 
this would require an accurate knowledge of the SFH from 
the redshift of the observations to the present. The uncer- 
tainty on the stellar ages and metallicities depends on the 
spectral signal-to- noise (S/N); a S/N ^ 20 ensures Alog(Z) 
< 0.3 dex, Alog(Age) < 0.2 dex and an average uncertainty 
on both parameters of ~0.12 dex. 

All galaxies in sample II which have a spectral S/N 
^ 20 and a stellar age and metallicity estimates computed 
as in lGallazzi et al.1 |2005l ) hereafter constit ute our sample 
S. Th is sample was extensively studied by iPasquali et al.l 
(|2010l ). who found that satellites are older and metal richer 
than centrals of the same M*, and this difference increases 
with decreasing M*. In addition, the average age and stel- 
lar metallicity of low-mass satellites (M* ^ 10 10 h~ 2 Mq) 
increase with the mass of the halo in which they reside. The 
authors interpreted these trends as due to the quenching of 
star formation in satellites, which leaves their stellar popu- 
lations to evolve passively (thus explaining their older ages), 
and to tidal stripping, which removes a non- negligible frac- 
tion of their stellar mass. In this picture, the stellar metal- 
licity of present-day satellites reflects the maximum stellar 
mass that they reached during their lifetime and which can 
be significantly higher than their present-day M*. 



2.2 Gas metallicities and star formation rates: 
sample G 

Sample II was also cross-matched with the catalogue of 
gas-phase meta l licitie s [i.e. 12 + log(0/H)] published by 
ITremonti et~ai1 l |2004 ) for the galaxie s in S PSS DR4. The 
analysis performed by ITremonti et all (|2004r i can briefly be 
summarized in the following steps: i) the continuum emis- 
sion in each galaxy spectrum was mod elled with a combi- 
nation of single stellar populations from lBruzual fc Chariot! 
(2003) of the same stellar metallic it y, also convolved with the 
extinction law of ICharlot fc F all (2000), and subsequently 
subtracted from the observed spectrum, it) The lines in 
the resulting pure emission-lines spectrum were simultane- 
ously fitted with Gaussians by imposing that all the Balmer 
lines have the same width and velocity offset, and likewise 
for the forbidden lines. Hi) The gas-phase metallicity was 
statistically estimated for all star forming galaxies by fit- 
ting the more prominent e mission lines with the models of 
ICharlot fc Long hctti (200l|) . These models combine th e syn- 
thetic stellar populations of iBruzual fc Charlotl ([2003) with 
the photoionization c ode CLOUDY (Ferland 1996) and the 
ICharlot fc Fall (|200(t) extinction law, so that the galaxy ISM 
conditions (e.g. gas ionization, metallicity and dust attenu- 
ation) are intertwined with the radiation field and nuclear 
synthe sis of the under l aying stellar populations. For each 
galaxy ITremonti et al. (2004) derived a likelihood distribu- 
tion of the gas-phase metallicity from which we extracted 
the median, the 16th and 84th percentile values. All galax- 
ies in sample II whose H/3 equivalent width is EW (H/3 ^ 
-3 A (corresponding to a S/N 10) and whose gas- phase 



metallicity is available from ITremonti et alj |2004l ) make up 
our sample G. 

For the galaxies in sample G we retrieved the median 
values of stellar mass, mflb, star formation rate, SFRfib, and 
specific SFRfib (sSFRflb = SFRfib/mfib) in the fibre, the me- 
dian values of global star formation rate, SFR g i , and spe- 
cific SFR glo (sSFR glo = SFRgio/M*) from the SDSS Data 
Release 7. For star forming galaxies (like those in sample 
G), the SFRflb was directly co mputed from their emissio n 
lines following the procedure of iBrinchmann et al. I (|2004h . 
who modelled the observed emission lines with the models of 
ICharlot fc Longhettil l|200ll ). The global star formation rate 
of the galaxies in DR7 was derived by fitting the photome- 
try of the outer region s of galaxies with th e synt hetic stellar 
population models of IBruzual fc Charlotl (12003) convolved 
with the extinction law of ICharlot fc Falll 1120001) (cf. Salim 
et al. 2007). Using this same procedure, Brinchmann et all 
l| 20041 ) derived the global and fibre stellar masses, M* and 
mfib, from the galaxy integrated magnitudes and the mag- 
nitudes measured within the fibre. 



2.3 Combined sample 

We extracted from sample G all galaxies for which the spec- 
tral S/N is larger than 20, and the stellar age and metallicity 
are available together with their gas-phase metallicity and 
star formation rate. We refer to this subset of objects as the 
composite sample, i.e. sample C. Table 1 summarizes the 
number of central and satellite galaxies in the samples S, 
G and C defined above. In Fig. 1 we compare the normal- 
ized distributions of the galaxies in sample II (black dashed 
line), sample S (grey solid line) and G (black solid line) in 
colour, concentration, stellar and halo masses. We note that 
the sample definition introduces a clear separation in colour 
and concentration between sample G and sample 5*. Being 
star forming, the galaxies in sample G turn out to be bluer 
and more disk-like (with lower conc entrations) than thos e 
in sample S. As already shown by IPasquali et al. I l|2009l ). 
star forming galaxies have on average lower M* and pref- 
erentially reside in low-mass haloes. This trend is mainly 
due to the central galaxies of sample G as it can be inferred 
from the right-hand panel of Fig. 2, while the satellites in 
the same sample span a wide range of environments from 
M h ~ 10 11 to 10 15 /i^Mq. In addition, centrals seem to 
be slightly more massive and metal-poorer than satellites in 
sample G (left-hand and middle panels of Fig. 2). It is also 
interesting to compare sample S and sample C in stellar 
age and metallicity; we find in Fig. 3 that the star forming 
galaxies in sample C have on average younger and metal- 
poorer stellar populations than those in sample S. We are 
thus probing different ranges of galaxy properties with re- 
spect to the analysis in Pasquali et al. (2010). 



3 THE STELLAR MASS - GAS-PHASE 
METALLICITY RELATION 

The dependence of gas-phase metallicity on stellar mass is 
plotted in Fig. 4 for centrals (grey) and satellites (black). 
We computed the median (solid lines), the 16th and 84th 
percentile values (dashed lines) of the gas-phase metallicity 
distribution of centrals and satellites in sample G per bin of 
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Figure 3. The normalized distributions of galaxies in sample S (grey) and sample C (black) in stellar age and metallicity. The distri- 
butions are not weighted by 1/V ma x. 



M^. It turns out that centrals and satellites follow a qual- 
itatively similar relation between gas-phase metallicity and 
M*, whereby the metallicity of their gas increases with their 
stellar mass. The gas in satellites is, though, metal richer 
than in centrals at nearly all stellar masses, and this is true 
for the median values as well as for the 16th and 84th per- 
centiles. 

The difference in the median 12 + log(0/H) be- 
tween satellites and centrals decreseas from 0.06 dex at 
log(M*/h- 2 M ) ~ 8.25 to 0.004 dex at log(M*/h" 2 M ) ~ 
10.75. This trend holds also for volume limited sub-samples 
extracted from sample G. 

For each bin of stellar mass, we constructed the cumula- 
tive distributions in 12 + log(0/H) of centrals and satellites, 
and used the Kolmogorov - Smirnov (KS), two sample test 
to assess the statistical significance of the metallicity offset 
seen in Fig. 4. The probability at which the null hypoth- 
esis (i.e. centrals and satellites are drawn from the same 
population) is rejected is larger than 95% in the range 8 
< log(M*/h~ 2 M ) < 10.5, and drops to ~19% for galax- 
ies more massive than M* = 10 10 ' 5 h~ 2 M0. Thus, over the 
stellar mass range spanned by our sample, the cumulative 
distributions indicate that satellites have a higher gas-phase 
metallicity than centrals. 

In order to check whether/how the M* - gas-phase 
metallicity relation of satellites depends on environment, 
we derived it for satellites split among different bins of 
halo mass. The results are shown in Fig. 5, where satel- 
lites residing in different groups are colour coded; the black 
dashed line traces the median 12 + log(0/H) of centrals 
from Fig. 4 and the grey band represents the 16th - 84th 
percentile range of the distribution of centrals also from 
Fig. 4. The la errorbars represent the error on the mean 
gas-phase metallicity weighted by 1/V m ax, and were com- 
puted by propagating the uncertainties associated with the 
individual measurements of 12 + log(0/H). Satellites less 
massive than rog(M»//i~ 2 M0) ~ 9.3 and in groups more 
massive than log(Mh/h~ 1 M.o) ~ 12 appear to follow the 
same M* - gas-phase metallicity relation. The overall shift 
to higher gas-phase metallicities of satellites with respect to 
centrals is due only to satellites in haloes more massive than 
10 12 /i _1 Mq, with their stellar mass - gas-phase metallicity 
relation slightly steepening as increases. 

Satellites in the least massive groups (11 < 
log(Mh//i _1 Me) < 12) and more massive than 
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Figure 4. The M* - gas-phase metallicity relation for central 
(grey) and satellite (black) galaxies. The solid lines represent the 
median metallicity, while the dashed lines the 16th and 84th per- 
centiles of the metallicity distribution of centrals and satellites in 
each bin of stellar mass. Centrals and satellites belong to sample 
G. 



log(M,// l - 2 M Q ) = 9.3 are very similar to centrals in 
terms of their gas-phase metallicity, while less massive 
satellites exhibit higher 12 + log(0/H) values than centrals 
as their stellar mass decreases [the maximum difference, 
+0.11 dex, is registered at log(M*/fr" 2 M ) ~ 8.2]. 



4 THE HALO MASS - GAS-PHASE 
METALLICITY RELATION 

We show in Fig. 6 the dependence of the gas-phase metallic- 
ity of central and satellite galaxies on the halo mass of the 
group in which they reside. The grey and black lines refer 
to centrals and satellites, respectively; the solid and dashed 
lines represent the 50th and 16th/84th percentiles of the 
gas- phase metallicity distribution within each bin of Mh- 
In halos with M h < 10 12 /i _1 M Q the gas-phase metallicity 
of central galaxies increases sharply with halo mass. Cen- 
trals in more massive halos, though, have roughly constant 
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Figure 5. The Af* - gas-phase metallicity relation for centrals 
(grey) and satellites (colour) split among different bins of halo 
mass. The black dashed line represents the median 12 + log(0/H) 
of centrals, while the grey band marks the 16th - 84th percentile 
range of the distribution of centrals. 



gas-phase metallicities. Note that the small number statis- 
tics for massive groups in the Yang et al. (2007) catalogue 
prevent us from probing this halo mass dependence beyond 
M h ~ 1O 14 /i _1 M . Regardless of halo mass, central galaxies 
are metal richer than satellites in halos of the same mass by 
~0.5 dex on average. This is simply a reflection of the fact 
that, by definition, central galaxies are the more massive 
galaxies in their groups, combined with the fact that gas- 
phase metallicity is an increasing function of stellar mass. 
The oxygen abundance of satellites smoothly increases with 
halo mass up to Mh — 1O 13 /i -1 M0, and is constant in more 
massive groups. Overall, it is seen increasing by ~0.5 dex 
across the full range of environm ents probed b y samp le G, 
in agreement with the findings of ICooper et a l. (2008). 

We split the satellites in sample G among different bins 
of stellar mass and computed their Mh - gas-phase metallic- 
ity relation as a function of M*. The results are illustrated 
in the left hand panel of Fig. 7, where satellites are colour 
coded on the basis of their stellar mass. The black dashed 
line traces the median 12 + log(0/H) of centrals as a func- 
tion of Mh (from Fig. 6) , while the grey band indicates their 
16th-to-84th percentile distribution in gas-phase metallicity 
as shown in Fig. 6. The la errorbars for satellites represent 
the error on the mean gas-phase metallicity weighted by 
1/V m ax, and were computed by propagating the uncertain- 
ties associated with the individual measurements of 12 + 
log(0/H). The offset in 12 + log(0/H) among the different 
Mh - gas-phase metallicity relations is simply due to stellar 
mass. We see that the gas-phase metallicity of satellites in- 
creases with their halo mass in all bins of M* . The amplitude 
of this trend across the full Mh range amounts to ~0.15 dex 
for satellites with log(M*//i~ 2 M Q ) < 9, ~0.1 dex for those 
with 9 < log(M*//i~ 2 M Q ) < 10.5 and only ~0.02 dex for 
the most massive satellites. These findings are consis tent 
with the results of lMouhcine. Baldrv fc Bamfordl (|2007l ). To 
better illustrate the dependence of gas-phase metallicity on 
M h , we computed the slope d[12 + log(0/H)]/d[log(M h )] in 
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Figure 6. The M h - gas-phase metallicity relation for central 
(grey) and satellite (black) galaxies. The solid lines represent the 
median metallicity, while the dashed lines the 16th and 84th per- 
centiles of the metallicity distribution of centrals and satellites in 
each bin of stellar mass. Centrals and satellites belong to sample 
G. 

each stellar mass bin, and plotted it as a function of M* in 
the right hand panel of Fig. 7. A weak trend can be recog- 
nized whereby the strength of the Mh - gas-phase metallicity 
relation decreases with increasing stellar mass. 



5 DISCUSSION 

In the previous sections we have shown that the gas-phase 
metallicity of central and satellite galaxies depends primar- 
ily on their ste ll ar ma ss. In agreement with the results of 
llremonti et alj (|2004l ). their gas becomes metal richer as 
their stellar mass increases. However, at fixed M* the gas- 
phase metallicity is higher in satellites than in centrals, and 
this observed offset increases with decreasing stellar mass. 
In addition, at fixed M* the oxygen abundance of satel- 
lites is found to be higher in more massive haloes. This 
trend is observed for satellites less massive than M* ~ 
lO lo ' 5 /i -2 M0 with an amplitude of ~ O.ldex over the range 
11 < \og(Mh/h~ M©) < 14, and becomes considerably 
weaker for the most massive satellites. Hence, there is a clear 
and 'pure' dependence of the gas-phase metallicity of satel- 
lites on environment, but this dependence is much weaker 
than the dependence on stellar mass. In what follows we in- 
vestigate the physical mechanisms that may be responsible 
for the (M*,Mh) - gas-phase metallicity relations of satellite 
galaxies. 

5.1 Gas-phase and stellar metallicities 

We start by comparing the properties in common between 
the galaxies in sample C and those in sample 5*. We re- 
call that sample C is built from sample S and includes only 
galaxies with gas-phase metallicity, stellar age and metallic- 
ity, while sample S contains galaxies with measured stellar 
age and metallicity, regardless of whether their gas-phase 
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Figure 7. Left: the Mj, - gas-phase metallicity relation for satellite galaxies (colour) split among different bins of stellar mass. The black 
dashed line represents the median metallicity of centrals, while the grey band shows the 16th - 84th percentiles distribution of centrals 
as a function of halo mass. Centrals and satellites belong to sample G. Right: the slope of the - gas-phase metallicity relations shown 
in the left hand panel is plotted as a function of stellar mass. The errorbars indicate the la uncertainty. 



metallicity is available (see Sect. 2). Similarly to sample G 
(cf. Fig. 4), also the satellites in sample C display a higher 
median oxygen abundance when compared with equally 
massive centrals. Figure 8 shows that their difference in me- 
dian 12 + log(0/H) at fixed M* decreases from ~0.14 dex 
at Mi, < 10 9 /i~ 2 M© to ~0.02 dex at Af* < 10 10 - 3 /i" 2 M Q . 
Note that this difference is larger than that measured be- 
tween centrals and satellites in sample G as shown in Fig. 
4. This is likely to be a consequence of the fact that sam- 
ple C is biased towards galaxies with weaker emission lines 
and lower star formation activity than those in sample G. In 
this respect, satellites in sample C may have been more sig- 
nificantly quenched by their environment than satellites in 
sample G. We caution, though, that sample C suffers from 
small number statistics when compared to sample G. The 
KS two sample test performed on the cumulative distribu- 
tions in 12 + log(0/H) of centrals and satellites in sample 
C in different bins of stellar mass establishes that satellites 
are metal richer than centrals at a confidence level between 
78% and 100% in the range 8 < log(Af*/ft- 2 M ) < 10.5, 
while the most massive satellites can not be distinguished 
from equally massive centrals. 

Using sample ff. lPasquali et all (|20 10T ) found that satel- 
lites have older, luminosity- weighted ages and higher stellar 
metallicities than equally massive centrals. Both these differ- 
ences become smaller for more massive galaxies. In addition, 
they showed that satellites with < 10 10 h~ 2 Mq have ages 
and metallicities that increase with the mas s of the host halo 
in which they reside. iPasquali et al.l (|2010l ) explained these 
trends by invoking i) strangulation, i.e. the removal of the 
gas reservoir of satellites after they are accreted by their 
host halo. The consequent quenching of star formation leaves 
the stellar populations of satellites to evolve passively, while 
the more prolonged star formation activity of centrals keeps 
their luminosity- weighted age younger, ii) Tidal stripping, 
which causes satellites to lose stellar mass once they have 
been accreted onto their host group. The fact that present 



9.2 



8.4 




_L 



Centrals 
Satellites 

ill 



_L 



9.5 10 
log(M,/h- 2 M ) 



10.5 



Figure 8. The M* - gas-phase metallicity relation for central 
(grey) and satellite (black) galaxies in sample C, i.e. with mea- 
sures of both gas and stellar metallicity. The solid lines represent 
the median metallicity, while the dashed lines the 16th and 84th 
percentiles of the metallicity distribution of centrals and satellites 
in each bin of stellar mass. 



day satellites descend from more massive progenitors is left 
imprinted in their stellar metallicity, which is not expected 
to change significantly after a satellite is accreted onto a 
bigger halo and its star formation activity is quenched. The 
effects of strangulation and tidal stripping likely depend on 
the time of infall, i.e. when a galaxy became a satellite for 
the first time. The earlier the infall time is, the older and 
less massive a satellite is, and the more massive the group in 
which the satellite resides (given that more massive haloes 
assemble earlier). This dependence on infall time naturally 
explains why the ages and stellar metallicities of satellites 
increase with M-^. 
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Figure 9. The distributions of stellar age and metallicity as a function of stellar mass for centrals (grey) and satellites (black) bcloging 
to sample S (left hand panels) and to sample C (right hand panels). The solid lines represent the median age/metallicity, while the 
dashed lines the 16th and 84th percentiles of the age/metallicity distribution of centrals and satellites in each bin of stellar mass. 



Since both stellar and gas-phase metallicities correlate 
with M*, one might also expect the higher gas-phase metal- 
licities of satellites (with respect to centrals of the same M*) 
to be a manifestation of tidal stripping. In this scenario, 
satellites have higher gas-phase metallicities than centrals 
of the same present-day mass, simply because they had a 
higher stellar mass at infall: the gas-phase metallicity serves 
as an indicator of the maximum stellar mass the galaxy 
reached during its lifetime. In order to test this hypoth- 
esis, we compare in Fig. 9 the M* - stellar age and the 
M t - stellar metallicity relations of centrals and satellites 
in sample S with those of centrals and satellites in sample 
C. As already pointed out by Fig. 3, the galaxies in sam- 
ple C are characterized by younger stellar ages and lower 
stellar metallicities than those in sample S. At fixed stel- 
lar mass satellites are systematically older than centrals in 
both samples. This age difference is rather small for sample 
C (about 0.04 dex at M* < lO lo /i~ 2 M ), and statistically 
significant only for galaxies more massive than lO 9 l5 ft- 2 M . 
In fact, the KS two sample test applied to the age cumula- 
tive distributions of sample C indicates that centrals and 
satellites with log(M*//i _2 M0) > 9.5 are drawn from two 
different populations at a >92% confidence level, while less 
massive centrals and satellites belong to the same parent 
population. In sample S the age difference between satel- 
lites and centrals is larger, and decreases from 0.19 dex at 
M* < 10 9 /T 2 M S to 0.02 dex at M* < lO lo - B /r 2 M . The 



KS two sample test applied to the age cumulative distribu- 
tions of centrals and satellites in different M* bins shows 
that the hypothesis of centrals and satellites belonging to 
the same population is rejected at a >99% confidence level. 
We may thus conclude that the satellites in sample C were 
accreted later and their star formation activity has been 
quenched less than the satellites in sample S. Consequently, 
they should have also experienced less tidal stripping, so 
that their stellar metallicities should be closer to those of 
equally massive centrals. This is indeed the case; the bot- 
tom panels of Fig. 9 show that satellites in sample C have 
stellar metallicities that are basically indistinguishable from 
those of centrals of the same stellar mass. The differences are 
less than 0.02 dex, implying that satellites in sample C have 
lost <10% of their pre-infall stellar mass. For comparison, 
the difference in the median \og(Z / Zq) between satellites 
and centrals in sample S decreases from ~0.3 dex at M* < 
10 9 /i~ 2 M e to <0.02 dex at M* > lO lo /i" 2 M , indicating 
that satellites in sample S have lost between ~75% (for a 
present day M* = 10 9 h~ 2 M e ) and ^20% (for a present day 
M* = 10 10 I5 /i~ 2 Mq) of their stellar mass at infall. To sum- 
marize, Fig. 9 indicates that satellites in sample C have not 
suffered significant tidal stripping (yet), but have undergone 
quenching of their star formation activity at some, small de- 
gree. Therefore, we conclude that stellar mass stripping can- 
not explain why satellites have higher gas-phase metallicities 
than centrals of the same stellar mass. 
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Figure 10. The dependence of the gas-phase mctallicity of cen- 
trals (grey) and satellites (black) in sample C on their stellar 
metallicity. The solid lines trace the median 12 + log(0/H), while 
the dashed lines the 16th and 84th percentiles of the gas-phase 
metallicity in each bin of stellar metallicity. 



In further support of this conclusion, we directly com- 
pare the gas-phase metallicities of centrals and satellites 
at fixed stellar metallicity in Fig. 10. If the same stellar 
stripping responsible for creating the (small) offset in stel- 
lar metallicity between star-forming centrals and satellites 
were the only responsible for the observed offset in gas- 
phase metallicities, we would expect the relation between 
gas-phase and stellar metallicity to be the same for central 
and satellite galaxies. Figure 10 shows that clearly this is 
not the case: at \og(Z / Zq) < -0.4 (corresponding to M* < 
1O 9,5 /i _2 M0, see Fig. 9) satellites have higher gas-phase 
metallicities than centrals of the same stellar metallicity. 
Rather, Fig. 10 suggests that the differences in gas-phase 
metallicity between centrals and satellites have their origin 
in something that is specific to the gas and that has a larger 
impact on less chemically evolved (hence preferentially low- 
mass) galaxies. The scatter between gas-phase and stellar 
metallicity has been shown to be associated with variations 
in stellar surface mass density, star formation rate and gas 
mass fraction (Gallazzi et al. 2005). In the next sections we 
test whether any of these parameters may indeed be asso- 
ciated with the difference in gas-phase metallicity between 
centrals and satellites. 



5.2 Surface mass densities 



iKauffmann et alj (2003) have shown that the star forma- 
tion activity of low mass galaxies depends more on their 
surface mass density than on their stellar mass. In a naive 
closed-box model for chemical evolution, this would imply 
that, at fixed M*, galaxies with a higher surface mass den- 
sity have transformed a larger amount of gas into stars and 
hence increased their gas-phase metallicity. For the galaxies 
in sample G we define the surface mass density in the fibre 
as Efib = mfib/ area ah, where areata is the physical area 
covered by the SDSS fibre (3" in diameter) at the galaxy 
distance. We use Efib as a proxy for the past star formation 



history that took place over the same (central) regions from 
which the gas-phase metallicities have been measured. The 
solid lines in the left hand panel of Fig. 11 plots the median 
Efib of centrals (grey) and satellites (black) as a function of 
their stellar mass; the dashed lines indicate the correspond- 
ing 16th and 84th percentiles. The difference in Efib between 
satellites and centrals is ~0.02 dex at most, and changes 
from being negative at log(Af*//i _2 M0) < 9.5 (Efib is lower 
in satellites) to positive at higher M*. In the right hand panel 
of Fig. 11 the median oxygen abundance (solid line) of cen- 
trals (grey) and satellites (black) is shown as a function of 
their surface mas s dens ity; as expected from the findings of 
IKauffmann et al l |2003h . the gas-phase metallicity increases 
with Efib for both centrals and satellites. When matched 
in Efib, satellites with log(Efib) < 8.5 display a higher gas- 
phase metallicity than centrals in the median and 16th and 
84th percentiles values. We therefore conclude that the off- 
sets in 12 + log(0/H) are not associated with any offsets 
in surface mass density, as expected in a simple closed box 
model if the different gas-phase metallicities of centrals and 
satellites owe to differences in their star formation histories. 



5.3 Gas mass fractions 

Similarly to Efib, the gas mass fraction can be considered 
an indirect measurement of the past star form ation activity 
under gone by galaxies. It has been shown by IZhang et al.l 
(2009) that the HI gas mass fraction of galaxies decreases 
along the gas-phase metallicity - stellar mass relation as both 
and the oxygen abundance increase. Moreover, at fixed 
stellar mass, gas - poor galaxies exhibit higher gas-phase 
metallicities. This poses the question whether the differ- 
ence in gas-phase metallicity between central and satellite 
galaxies reflects a difference in their gas mass fractions. In 
order to test this, we computed the gas mass fraction in 
the fibre for central a nd satellite galaxies using Eq. (5) of 
iTremonti et alj <|2004T l ; i.e., we derive the surface gas mass 
density, E gas , from the star formation surface density in 
the fib re by inverting the Schmidt law given by iKennicuttl 
(1998). The gas mass fraction in the fibre is then defined 
as /Xgas = Egas /(Egas + Efib), where Efib is the surface mass 
density calculated in Sect. 5.2. The gas mass fraction is plot- 
ted as a function of Mi, in the left hand panel of Fig. 12, 
where centrals are in grey and satellites in black. Solid lines 
represent the median n sas in bins of stellar mass, while the 
dashed lines indicate the the 16th and 84th percentiles. The 
gas mass fraction decreases as galaxies become progressively 
more massive, and no significant difference is seen between 
centrals and satellites. We also determined the 16th, 50th 
and 84th percentiles of the gas-phase metallicity of centrals 
and satellites in different bins of ^ gas , as shown in the right 
hand panel of Fig. 12. For both centrals and satellites, the 
oxygen abundance decreases with incr easing gas mass frac - 
tion, in agreeme nt with the findings of IZhang et al.l (2009), 
lErb et alj ||2006l ) and Erb et al. (2008, for galaxies at z ~ 
2). However, at fixed ^i g a S satellites are typically metal richer 
than centrals, by 0.06 dex at /i gas — 0.1 to 0.12 dex at ^ gas — 
0.6. This trend likely arises from the anti-correlation be- 
tween ^tga S and Mi,, combined with the fact that, at fixed 
Mi,, satellites are metal richer than centrals. We conclude 
that the difference in gas phase metallicity between centrals 
and satellites is not associated with a difference in (indirectly 
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Figure 11. Left: the surface mass density in the fibre as a function of stellar mass for centrals (grey) and satellites (black). The solid 
lines show the median Egb while the dashed lines the 16th and 84th percentiles in different bins of stellar mass. Right: the dependence 
of gas-phase metallicity on £fl b for central (grey) and satellite (black) galaxies. 
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Figure 12. Left: the gas mass fraction in the fibre as a function of stellar mass for centrals (grey) and satellites (black). The solid lines 
show the median /i ga s while the dashed lines the 16th and 84th percentiles in different bins of stellar mass. Right: the dependence of 
gas-phase metallicity on /%as for central (grey) and satellite (black) galaxies. 



inferred) gas mass fraction, as would have been expected in 
a simple closed box model. 

5.4 Star formation rates 

The results in the previous two subsections show that the 
enhanced gas-phase metallicities of satellite galaxies are not 
associated with any enhancement or suppression of the stel- 
lar surface mass density or the gas mass fractions. We now 
test whether there is any correlation with the star formation 
rates. 

Fig. 13 shows the specific star formation rate in the fi- 
bre [sSFRflb = log(SFRflb/mfib)] and the specific global star 
formation rate [sSFR g i D = log(SFR g i /A/*)] as a function 
of stellar mass for centrals (grey) and satellites (black) in 
sample G. The sSFRflb of satellites is higher by only 0.02 
dex than that of equally massive centrals in the range 9 < 
log(M*//i~ 2 M©) < 10.5. The KS two sample test applied 
to the cumulative distributions in sSFRflb of centrals and 



satellites in different bins of stellar mass indicates that such 
a difference is not statistically significant Q We checked that 
the differences in gas-phase metallicity between satellites 
and centrals do not disappear when galaxies are matched 
in sSFRflb and stellar mass. 

Interestingly, the median sSFR g i of satellites is system- 
atically lower than that of centrals by 0.03 dex on average 
at any stellar mass. The KS two sample test run on the 
cumulative distributions in sSFR g i of centrals and satel- 
lites in different bins in M* indicates that this difference is 
statistically robust at >95% confidence level, at any stellar 
mass. The fact that satellites have lower global star forma- 
tion rates than centrals of the same stellar mass is consistent 
with numerous other studies (e.g., Weinmann et al. 2006, 
2009; Kimm et al. 2007; Cooper et al. 2007; van den Bosch 
et al. 2008a,b; Pasquali et al. 2010; Peng et al. 2010; Wetzel, 

1 The same result is obtained when central and satellite galaxies 
are compared in their sSFRflb defined as log(SFRflb/Af*). 
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Figure 13. The dependence of the specific star formation rate in the fibre (sSFRgb) and the specific global star formation rate (sSFR g i ) 
as a function of stellar mass for centrals (grey) and satellites (black) in sample G. The solid lines indicate the median specific SFR and 
the dashed lines the 16th and 84th percentiles of the distribution in specific SFR in different bins of stellar mass. 



Tinker & Conroy 2011), and indicates that certain processes 
operating on satellite galaxies causes them to quench their 
star formation. The fact that our results indicate that satel- 
lites have lower sSFR g i than centrals, but are equally active 
in their central regions (within the fibre) , suggests that this 
quenching is an outside-in process. 

5.5 Beyond a closed box 

The results in the previous subsections have shown that the 
differences in gas-phase metallicities between centrals and 
satellites cannot be understood within a simple closed-box 
model for their chemical evolution. This should not come 
as a surprise, since galaxies are not expected to be closed 
boxes. Rather they experience inflows, outflows and even 
mass stripping. Interestingly, the efficiency of these processes 
is expected to depend strongly on environment, which, as 
we argue below, has the potential of explaining the trends 
identified in this paper. 

5. 5. 1 Strangulation 

Central galaxies that reside in host haloes with 
M h < 1O 12 /i _1 M0 are expected to accrete cold gas 
from their surroundings ('cold flows'), whereas centrals in 
more massive haloes are expected to accrete most of their 
gas through cooling flows from their hot gaseous haloes 
(e.g., Birnboim & Dekel 2003; Keres et al. 2005). Satellite 
galaxies, however, are likely to be 'deprived' of these inflows, 
simply because their gas reservoir, be it hot or cold, is likely 
to be stripped by tidal forces and/or ram-pressure from the 
hot gaseous atmosphere associated with the host halo. Since 
inflows are likely to be metal-poor, this 'strangulation' 
(Larson et al. 1980) of satellite galaxies is expected to 
suppress the dilution of their ISM, resulting in gas-phase 
metallicities that are higher than those of central galaxies, 
in qualitative agreement with the data. Furthermore, as 
already discussed in Pasquali et al. (2010), present-day 
satellites in more massive haloes were, on average, accreted 
earlier (see also De Lucia et al. 2012). Hence, satellites in 
more massive haloes have been 'deprived of dilution' for a 



longer time than satellites of the same stellar mass in less 
massive haloes, which might explain why the gas-phase 
metallicity of satellites increases with host halo mass (see 
Fig. 7), at least qualitatively. 

However, strangulation is unlikely to be the entire pic- 
ture. First of all, if the differences in gas-phase metallicity 
(within the fibre) are due to differences in the amount of 
dilution, this should probably show up as differences in the 
gas mass fraction (again within the fibre). As is evident from 
Fig. 12, such difference is absent, although we caution that 
these 'gas mass fractions' are indirectly inferred from the 
star formation rates (within the fibre). We conclude that 
strangulation is likely to play a role in regulating the gas- 
phase metallicities of satellite galaxies, but that direct data 
on the gas mass fractions of satellite galaxies, combined with 
detailed hydrodynamical simulations of strangulation, are 
needed to test whether strangulation in itself can explain 
the various trends presented in this paper. 

5.5.2 Ram-Pressure Stripping and Galaxy Harassment 

In addition to strangulation, there are two other mecha- 
nisms that may impact gas-phase metallicites and which 
operate on satellites but not on centrals: ram-pressure strip- 
ping (e.g., Gunn & Gott 1972) and galaxy harassment (e.g., 
Farouki & Shapiro 1981; Moore et al. 1996). We now dis- 
cuss how these two processes may cause satellite galaxies to 
have higher gas-phase metallicities than centrals of the same 
stellar mass. 

An important mechanism for regulating the gas-phase 
metallicities of (the inner regions of) disk galaxies is ra- 
dial gas flows. Angular momentum redistribution within the 
disk, due to torques from bars, spirals and/or interacting 
galaxies can cause the cold gas at large galactocentric radii 
to move inwards (cf. Roskar et al. 2008; Minchev et al. 2012). 
Since galaxies typically have metallicity gradients that de- 
cline outwards, such radial gas flows will dilute the gas in the 
inner regions. Since satellite galaxies are exposed to stronger 
tidal forces, and to high-speed impulsive encounters (i.e., 
galaxy harassment), they are more likely to experience such 
diluting radial flows than central galaxies. If this were an 
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important mechanism, satellite galaxies would have central 
gas-phase metallicities that are lower than central galaxies, 
opposite to what is observed. This suggests that, contrary to 
what we postulated above, satellite galaxies do not experi- 
ence (an enhancement in) radial gas flows. This could come 
about if ram-pressure and/or tidal forces strip the outer, 
metal-poor gas before it has the opportunity to migrate in- 
wards. In that case, satellite galaxies are less likely to expe- 
rience diluting, radial gas flows than centrals, which again 
might explain why the latter have lower gas-phase metallici- 
ties th an t he former (see also Shields. Skillman fc Kennicuttl 
i|l99ll ) and lPetropoulou et al] (|201ll )'). 

There is ample observational evidence to support the 
notion that ram-pressure stripping is an important mech- 
anism, at least in cluster environments (e.g., Gavazzi et 
al. 2001; Solanes et al. 2001; Koopmann & Kenney 2004; 
Abramson et al. 2011). This is also supported by detailed hy- 
drodynamical simulations (e.g., Abadi et al. 1999; Roediger 
& Hensler 2005; Kronberger et al. 2008; Bekki 2009). How- 
ever, the overall efficiency of ram-pressure stripping (i.e., 
what fraction of the gas is stripped and on what time scale) 
is still a matter of debate, mainly because the simulations 
have shown that it is a strong function of the porosity of 
the galaxy's ISM (e.g., Quilis et al. 2000; Tonnesen & Bryan 
2009). What is well established, though, is that ram-pressure 
stripping is more efficient for less massive galaxies in more 
massive environments (e.g., Gunn & Gott 1972; Bekki 2009). 
This makes the simple prediction that the suppression of 
metallicity-dilution in satellite galaxies increases with de- 
creasing stellar mass and with increasing host halo mass; 
this is consistent with the observed trends, at least quali- 
tatively. Another effect that may contribute here is again 
the distribution of infall times; since low mass satellites are 
typically accreted earlier than their massive counterparts 
(e.g., De Lucia et al. 2012), they have been exposed to the 
'corrosion' exerted by their environments for a longer pe- 
riod. If ram-pressure stripping indeed causes satellite galax- 
ies to have gas-phase metallicities that are higher (i.e., less 
diluted) than those of centrals of the same stellar mass, then 
the statistics of infall times would contribute to explaining 
the observed trends of gas-phase metallicities of satellites as 
a function of M* and M^. 

An interesting consequence of ram-pressure stripping, 
which may have important consequences for our discussion, 
is that the gas that is not stripped away, which is mainly 
the gas at small galactocentric radii, may actually be com- 
pressed by the ram-pressure, giving rise to (significantly) 
enhanced star formation in the disk (e.g., Dressier & Gunn 
1983; Gavazzi et al. 1995; Kronberger et al. 2008; Kapferer 
et al. 2009). Our results that satellite galaxies have the same 
SFRs in their central regions (within the fibre) as centrals 
of the same stellar mass (see Fig. 13), therefore seems to 
argue that ram-pressure stripping is not very efficient. On 
the other hand, there are clear indications that the global 
SFRs of satellites are suppressed with respect to those of 
centrals, which might be an indication that ram-pressure 
stripping has started to remove gas from the outer disks. A 
ram-pressure stripping induced boost in the SFR will also 
lead to an increase in the amount of metals returned to 
the ISM. However, whether these metals can be retained or 
whether they are ejected and/or stripped is likely to depend 
on many factors, to the extent that it is unclear whether 



this boost in SFR is expected to increase or decrease the 
gas-phase metallicities. More detailed studies, in particular 
of the HI content of satellites versus centrals, will be needed 
to investigate the potential impact of ram-pressure stripping 
for explaining the gas-phase metallicities of satellite galax- 
ies. 



5.5.3 Galactic Wind Confinement 

Another environment-dependent process that may poten- 
tially play a role in 'regulating' the gas phase metallicities 
of satellite galaxies is galactic wind confinement. When a 
galaxy, central or satellite, is located within a dark mat- 
ter halo that has a hot gaseous atmosphere, the associated 
pressure may be able to inhibit a potential supernova-driven 
galactic wind from escaping the galaxy. Since galactic winds 
are typically metal-enhanced (e.g., Mac- Low & Ferrara 1999; 
Dalcanton 2007) , systems in which the winds can escape are 
expected to have lower gas-phase metallicities than systems 
in which the wind material is confined to be recycled within 
the ISM of the galaxy itself. Since satellite galaxies reside in 
more massive host haloes than central galaxies of the same 
stellar mass, galactic wind confinement might explain why 
the former have higher oxygen abundances than the latter. 

This mechanism also nicely explains the observed trends 
with both stellar mass and halo mass, at least qualitatively. 
After all, the hot gaseous atmosphere in more massive haloes 
has higher pressure, and is therefore better able to con- 
fine winds. This might explain why the difference in gas- 
phase metallicities between centrals and satellites is larger 
for satellites in more massive host haloes. In addition, since 
winds are expected to be more efficient in less massive 
haloes, which host less massive centrals, the potential im- 
pact of winds is more important for less massive galaxies, 
providing a qualitative explanation for the scalings with stel- 
lar mass. The recent simulations developped by Bahe et al. 
(2012) indicate that ram pressure is larger than the con- 
finement pressure independently of halo mass. Only a small 
fraction of the simulated galaxies (~16%) appears to be con- 
finement dominated, but this occurs when they have already 
gone through the first pericentre passage, where ram pres- 
sure is highest, and have hence lost most of their hot gas. 



6 CONCLUSIONS 



We used the groups catalogue o f lYang et all (2007) together 
wi th the DR4 catalog ue of stellar ages and metallicities 
of iGallazzi et al] (12005), the DR4 com pilation of gas- phase 
metallicities of Tremonti et al] ( 2004) and the DR7 cata - 
logue of star formation rates of iBrinchmann et al. I l|2004h . 
to study the gas-phase metallicity of central and satellite 
galaxies and its dependence on their stellar mass and the 
dark matter mass of their host environments. We find that: 

i) satellites have on average higher gas- phase metal- 
licities than central galaxies of the same M*. The magni- 
tude of this differences depends on stellar mass, increasing 
from 0.004 dex at log(M*//i" 2 M ) ~ 10.75 to 0.06 dex at 
log(M*//i- 2 M Q ) ~ 8.25. 

ii) Within the same host halo, satellites have lower gas- 
phase metallicities than their central galaxies by ~0.5 dex 
on average. This simply reflects that, by definition, central 
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galaxies are the most massive galaxies in their halos, com- 
bined wih the fact that gas-phase metallicity increases with 
stellar mass. 

iii) At fixed M*, the gas-phase metallicity of satel- 
lites increases with halo mass; this trend is more pro- 
nounced for less massive galaxies. Low mass satellites with 
W s - 5 h~ 2 M Q ^ Mi, ^ 1O 9 /i" 2 M have average oxygen 
abundances that increase by ~0.15 dex between Mh = 
1O 11 /i- 1 M and M h = lO 14 ft _1 M . 

iv) Satellite galaxies have higher gas-phase metallici- 
ties than central galaxies with the same stellar metallicities. 
The magnitude of this difference increases with decreasing 
metallicity. 

We have contrasted the gas-phase metallicities of cen- 
trals and satellites with other galaxy properties, such as stel- 
lar age, surface mass density, gas mass fraction and specific 
star formation rate, in order to identify the mechanism(s) 
responsible for the observational results listed above. Based 
on these comparisons we draw the following conclusions: 

• Star-forming centrals and satellites have very similar 
stellar metallicities at fixed M*. This basically rules out 
that the difference in gas-phase metallicities between cen- 
trals and satellites of the same stellar mass is a consequence 
of stellar mass stripping. Interestingly, this mechanism was 
used by Pasquali et al (2010) in order to explain a simi- 
lar difference between centrals and satellites but in stellar 
metallicity for mostly quiescent galaxies. The fact that this 
same mechanism cannot explain the gas-phase metallicities 
is also immediately evident from the fact that satellites have 
higher gas-phase metallicities than centrals of the same stel- 
lar metallicity (see Fig. 10). Rather, these results suggest 
that the differences in 12 + log(0/H) have their origin in 
something that is specific to the gas. 

• Star forming satellites and centrals show very similar 
surface mass densities and gas mass fractions (both mea- 
sured within the fibre) at fixed M*. Consequently, their ob- 
served offset in gas-phase metallicity remains the same when 
they are matched in surface mass density or gas mass frac- 
tion, rather than stellar mass. Since these quantities are an 
indicator of the galaxy's past, integrated star formation his- 
tory (SFH) within the fibre we conclude that the differences 
in 12 + log(0/H) between centrals and satellites (which are 
also measured within the fibre) are not a consequence of cen- 
trals and satellites having experienced different closed-box 
star formation histories. Rather, we argue that the mecha- 
nism^) responsible for this offset must be associated with 
differences in their inflow and/or outflow histories. 

• Despite their difference in gas-phase metallicities, cen- 
trals and satellites of the same stellar mass have present-day 
star formation rates, measured within the fibre, that are in- 
distinguishable. However, when corrections are made to in- 
clude the star formation occuring on larger scales, outside 
the fibre, we find that satellites have slightly reduced global 
SFRs compared to centrals of the same M*. In addition, we 
find that star forming satellite galaxies are characterized by 
somewhat older stellar ages than equally massive star form- 
ing centrals. These results suggest that the satellite galaxies 
have started an outside-in process of star formation quench- 
ing. We emphasize that the overall impact of quenching is 
small in our sample, which is a consequence of the fact that 
our sample is biased towards star forming galaxies, for which 



gas-phase metallicities can be measured from their emission 
lines. 

We argue that there are three mechanisms that oper- 
ate on satellite galaxies and that can potentially explain 
why they have gas-phase metallicities that are higher than 
those of central galaxies of the same stellar mass: (i) stran- 
gulation, which inhibits the inflow of metal-poor gas, (ii) 
ram-pressure stripping of the outer, metal-poor gas-disk, 
thereby preventing dilution of the central gas due to radial 
gas flows promoted by torques, and (iii) pressure confine- 
ment of metal-enriched outflows. As we have discussed in 
the text, each of these three mechanisms naturally explains 
why the difference in gas-phase metallicity between centrals 
and satellites increases with decreasing stellar mass and with 
increasing host halo mass, at least qualitatively. More work 
needs to be done in order to quantitatively discriminate be- 
tween these environmental processes. For example, detailed 
simulations of ram pressure stripping and galactic wind con- 
finement, which include also star formation and stellar/AGN 
feedback, are required in order to predict the effect of these 
mechanisms on the observed stellar and gas-phase metal- 
licities of satellites. From an observational point of view, 
determining the star formation history of satellites of differ- 
ent stellar mass and in different environments will constrain 
the time-scale over which strangulation/ram pressure strip- 
ping quench star formation, while a direct measurement of 
the (central and total) gas mass of satellites will test the 
efficiency of strangulation and ram pressure stripping. 
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